function [der0]=read_data(data_name,n_nodes)

file_name=strcat(data_name,'.dat_line');
tem=fopen(file_name,'r');
der0=fread(tem,n_nodes,'real*8');
fclose(tem);